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Abstract 

We consider the rare decay t] — > 7r°7r°77 and calculate the non-resonant contribution 
to the amplitude to one loop in Chiral Perturbation Theory. We display our result 
as both a diphoton energy spectrum and a partial decay rate as a function of the 
photon energy cut. It turns out that the one-loop correction can be numerically very 
important and could be detected, at sufficiently large center-of-mass photon energies, 
from a measurement of the partial decay width. 
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1 Introduction 



The 77 — > ^7777 decays been been recently analyzed in the framework of Chiral 
Perturbation Theory (CHPT) [0, H (see 0] for a recent review). The experimental 
interest for such rare decays stems from the large number of observed 77's anticipated 
at various //-factories, e.g. CELSIUS, ITEP and DA$NE ||, as well as at other 
facilities, such as GRAAL, MAMI, ELSA and CEBAF. 

In the neutral decay 77 — > 7r°7r°77 we can identify two physically distinct contri- 
butions to the amplitude: 

A(t] -> TrVV) = A R + A NR . (1 

The first contribution (that we call the 'resonant amplitude') is characterized by the 
7T° pole in the diphoton invariant mass squared (s 77 ) and is proportional to the on- 
shell amplitudes 77 — > 37r° and 7r° — > 77:^ 

A{r] _> 3n°)A(n° -» 77) n 

^ifi = — — 2 • y 1 - 

By construction, Ar can be predicted up to a phase from the experimental data on 
77 — ► 37T° and 7r° — > 77 and needs not to be calculated in CHPT. The main result of PJ, 
in the 77 —>■ 7r°7T°77 channel, can be summarized as follows: the resonant contribution 
dominates over the tree-level non-resonant one in the whole phase space. 

The non-resonant amplitude cannot be predicted using experimental data and 
must be calculated in CHPT. In Ref. [|]] Anr has been calculated at the tree level, 
i.e. to 0(p 4 ); at this order only the 77-exchange diagram (77 — > 7t tt 77* — > 7r 7r°77) 
contributes. Knochlein et al. have considered also the //-exchange diagram (formally 
of higher order in CHPT), but for both 77- and //-exchange they found that A^r is 
negligible with respect to Ar. This is because the lowest-order 7?7?7r 7r and rjrj' Tr°Tr° 
vertices vanish in the limit m u = = 0. In spite of the analogous suppression factor 
in the 7r°-exchange contribution, which is proportional to m u — ma, the enhance- 
ment due to the pole makes Ar dominant with respect to A^r over the full range of 
kinematical parameters. However, there are good reasons to presume that this sup- 
pression does not occur to one loop. Indeed, recently Talavera et al. have shown 
that the one-loop contribution to the 77 — ► 7r 7r°7r° amplitude - which is obviously 
connected to the 77 — > 7r°7r°77 one - dominates over the tree-level result (the former 
is one order of magnitude larger than the latter). The reason for this enhancement 
can be traced to the fact that the tree-level amplitude is proportional to m^, whereas 
the one-loop correction is not affected by such suppression. 

In this Letter we compute A NR to one loop, neglecting isospin-breaking effects 
and the suppressed 77-exchange diagrams. In this limit only one-particle-irreducible 
(1PI) diagrams contribute and the result is finite. Furthermore, A^r does not receive 
any contribution from the 0(p 6 ) counterterms analyzed in 0, just as it happens in the 
77 — > 7r°7r°7r case. We find that the one-loop contribution to A^r is dominant, with 



1 The minus sign in Eq. (1.2) is due to our convention for the amplitudes. 
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respect to the corresponding tree-level one. Once again, the reason is that the tree- 
level amplitude goes to zero in the limit m u = = 0, whereas the one-loop result 
does not vanish in this limit. We find, in addition, that there exists a kinematical 
region in the phase space, i.e. a region of sufficiently large values of s 77 , where the 
resonant amplitude - which is a background that shadows A^r - is suppressed, and 
a measurement of the decay width T(r] — > 7r°7r°77) would allow us to detect a purely 
0(p 6 ) effect. 

The outline of this Letter is as follows. We begin in Section 2 with the description 
of the kinematical variables for the decay rj — > 7r°7r°77 and a list of interaction terms 
relevant to the calculation of the one-loop amplitude. In Section 3, after recalling 
the expression of the decay amplitude at the tree level, we describe the calculation 
of the one-loop corrections to the non-resonant amplitude and give the analytic 
expression of the result. We proceed in Section 4 to calculate the decay width, starting 
from the sum of A R + A NR . We display the result in the form of both a diphoton 
energy spectrum and a partial decay rate as a function of the energy cut around 
si,/ — m w ■ Then we determine the phase-space region, in terms of ere, where 

the suppression of the background due to A R may allow to detect the one-loop effects 
in A NR . We end the Letter with some concluding remarks and a short discussion of 
further developments, including possible extensions of the one-loop calculation to the 
charged pions channel, as well as to the reactions 77 — > r/nn. 



2 Kinematical variables and interaction terms 

The kinematics of the decay r](q) — > 7r°(pi)7r° (2)2)7 ^1)7(^2, £2) can De described in 
terms of five independent scalar variables which we choose as: 

Stttt = (Pi + P2) 2 , 21,2 = h,2 ■ (Pi + P2) , (2.3) 

s 77 = (£4 + k 2 f , z 3 = (£4 + k 2 ) ■ (pi - p 2 ) ■ (2.4) 
If we write the decay amplitude in the following way 

A( v - 7r°7r°77) = e 2 ^A^ , (2.5) 

the decay width is given by 



Since the process i] — > 7r°7r°77 involves the electromagnetic interaction of an odd 
number of pions, the decay amplitude receives contributions only from the odd- 
intrinsic parity sector of CHPT and thus is at least 0(p A ). 

The CHPT lagrangian, expanded up to 0(p 4 ), is given by 

C = £ (2) + £ (4) , (2.7) 
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where 

F 2 

£( 2 ) = _tr (D^UD»U ] + X U ] + X ] U) (2.8) 

and can be split into the odd-intrinsic anomalous part (i.e. the Wess-Zumino 
term fH) and the 0(p 4 ) Gasser-Leutwyler lagrangian || 

10 

= C wz + £ ULf . (2.9) 

i=i 

As usual, we assume the exponential parametrization U = exp (iV2P 8 /F), where P 8 
is the SU (3) octet matrix of pseudoscalar mesons and F coincides to the lowest order 
with the charged pion decay constant F n = 92.4 MeV |J. The covariant derivative 
in Eq. (|2.8|) is given by D^U = d^U + ieA^[Q, U], where is the photon field and 
Q =diag(2/3, —1/3, —1/3). Finally, we employ the identification x = = 2BA4 in 
the external scalar sources, where M. = diag(m u , rrid, m s ) is the quark mass matrix 
and B can be identified to the lowest order with the mass ratio B = m 2 /(m u + rrid). 

In principle, in this decay one could take into account also the 77 — 77 mixing, i.e. 
the mixing of Pg with the singlet-field 1] . However, as shown in Jl|, this effect can 
be safely neglected in the non-resonant amplitude. Hence in the loop calculation we 
identify the mass-eigenstate 77 with the octet field 7?g. In the same spirit, since we 
neglect isospin-breaking effects in Anr, in the following we assume m n = m^o. 

For the interaction terms necessary to calculate the tree-level amplitudes one can 
refer to [|TJ. However, not just for the sake of completeness, but also given that we 
use a small subset of the couplings in D, we collect them in the following: 

A^( m ^nW) = 3AW(r ]8 ^n\+n-) = Bo{m "- md) , (2.10) 
A^s-^V) = A^\r ]s ^r ]8 n^-) = Boim 3 u F + 2 md) , (2.11) 

A (4) (7r o^ 77) = v^ w (ift^Ty) = 4^r^ef^^. (2.12) 

In addition, in order to compute the one-loop diagrams in fig. [I], we introduce the 
generic couplings 

A m (0+0- ^ 0000) = as ^ + fom 2 + ^ _ m 2) + d{p 2_ _ m 2 } ^ (2 13) 

A (4) ( 0O^0+0- 7) = fe^Wjjtf , (2.14) 

where q, p±, pi )2 and k are the (outgoing) momenta of the pseudoscalars 0°, ± , ^2 
and of the photon, respectively. The constants a, b, c and d have the dimensions of 
inverse mass squared, whereas / has those of an inverse mass cubed. As we will show 
in the next section, the 'off-shell couplings' c and d are irrelevant for our process, 
since their contribution to the amplitude cancels out as a consequence of the gauge 
invariance. In the tt + tt~ — > 7r 7r° and r/ 8 — > 7r + 7r~7 cases, useful in order to estimate 
the dominant pion loops, we find 

and /= -i^- (2 - 15) 
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Figure 1: 1PI one-loop diagrams for the (jr — > 0^0277 transition. The second photon line 
has to be attached to the charged lines running in the loop and to the vertices. 

3 The decay amplitude to one loop 

In this Section we would like first to briefly recall the expression of the tree-level 
amplitudes Ar and A^ R , obtained by considering the 7r° and r] 8 exchange diagrams 
1 

4 (4) _ e 2 B (m u -m d ) ^ h u p « h P ( o-i\ 



4V3?r 2 F3 (s 



,(4) _ e 2 E (m M + m d ) e n k v e a k P ( o 9) 

NR ~ YZV^^^W" 

As already stated in the introduction, the enhancement factor due to the pion pole 
makes A$ dominant -in spite of its suppression factor (m u — m d )- with respect to 
Atf R in the entire kinematical space. 

The 0(p 6 ) loop and counterterm (CT) contributions can be divided in three 
gauge-invariant subgroups: reducible 7r°-exchange diagrams, reducible r/8 _exc hange 
diagrams and 1PI diagrams. 

i. The 7r°-exchange diagrams, which include both loops and CT, contribute mainly to 
Ar. In principle such diagrams generate also a contribution to A^r. Indeed, if we 
decompose the r] — > 7r°7r°(7r )* amplitude as follows: 

A(r] -> 7r°7r°(7r )*) = A on _ shell (7] -> 3tt°) + (s 77 - m 2 ) x A off _ shell , (3.3) 

the A ff -shell term drops out of Ar. However this non-resonant contribution vanishes 
in the limit m u = rrid and thus can be safely neglected. On the other hand, \A on -sheii\ 
can be extracted by means of experimental data and we do not need to evaluate it in 
CHPT. 

ii. Also the ^g _exc hange diagrams include both loops and CT. These diagrams contribute 
only to A NR and can be safely neglected. Indeed, we have explicitly checked that the 



contribution of these diagrams is of the same order as the tree-level result ( |3.2| ) that is 
known to be small Jl| . The reason of such suppression can be easily understood: the 
7T-7T loops, which are expected to provide the dominant contribution, are suppressed by 
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the factor (m u + md) in (|2.11|) just as the tree level result. There are contributions from 
the K-K loops and which are not suppressed by (m u + m^) . They are nonetheless 
negligible, since we are far below the kaon threshold and the CT combinations involved, 
i.e. (Li + L 3 /6), {L 2 + L 3 /3) and L 4 , are small ||. 

iii. The 1PI diagrams are the loop diagrams in fig. 1 (note that the figure indicates schemat- 
ically at least four distinct diagrams). The sum of these contributions is finite and, as 
we will show in the following, turns out to be the dominant contribution to Anr. 

The calculation of the loop diagrams in fig. 1 resembles a recent calculation by 
D'Ambrosio et al. of the radiative four-meson amplitudes. The difference is that 
in our case one pseudoscalar field is replaced by one photon, but the main features of 
the result -simply dictated by QED- are the same. We find indeed 



AlPI 

^NR 



4ef{as nn + bm\ 



x <j C 20 (s n7T , -z 2 )e^ a pe i lk u l [(e 2 • pi 2 )k 2 - <f + ( e b fa «-> e 2 , fa) \ , (3.4) 



where p\ 2 = Pi + P2 and the function C 20 (x,y) is given in appendix. As in (10 
gauge invariance forces the result to depend only from the 'on-shell couplings' a, b 
and /. Furthermore, the amplitude (|3.4j) is 0(fa, fa) in the limit of vanishing photon 



momenta, in analogy to the direct-emission amplitudes of |TD[ which are 0(k). 

Since we have written the two vertices in a general form, not only the dominant 
pion loops, but also the kaon loops are represented in Eq. fl3.4|) . From our general 
one-loop amplitude we recover, as a particular case, part of the result of Talavera et 
al. H (i.e. the contribution of 1PI diagrams). The precise correspondence between 
the function C 20 (x,y) and the function R(x,y) entering Eq. (10) of is given by 



R{x,y) = 32ir 2 yC 20 (x,y) 



(3.5) 



The amplitude ( |3.4j) depends only on the function C 2 $ and thus is finite. This is 
a consequence of both the gauge invariance of the amplitude and the fact that the 



on-shell ti + tt 



7r°7r° 



amplitude depends only on (i.e. it does not depend on the 



loop variables). We expect that the sum of 1PI diagrams is no more finite if the two 

pair. Indeed, in this case not only the on-shell 



7T + 7r amplitude does depend on the loop momenta, but this sum is also not 



external tt s are replaced by a 7r -7r 

7T + 7T~ 

gauge invariant (in order to obtain a gauge invariant result, it is necessary to add the 
corresponding reducible diagrams with a photon emission from the external legs). 



4 Numerical analysis 



The results of our analysis are summarized in figs. 2 and 3. The plots have been 
obtained integrating numerically Eq. (|2.(j|) with the following decay amplitude: 



A(r] -> 7r°7r°77) = A p R hys + 



4(4) , AlPI 

NR ^ ^NR 



(4.1) 
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Figure 2: Diphoton spectrum (zy = s 77 /m^) for the decay 77 — > 7r°7r°77. The upper full 
line is the total contribution. The dashed line is the resonant contribution (lyr^^l 2 ), the 
dotted line is the one-loop non-resonant contribution (|A^| 2 ) and the dash-dotted line is 
the asbolute value of their interference (p = 2, «o = 0.18). The lower full line is the tree-level 
non-resonant contribution (|t4^| 2 ). 



The last two terms in Eq. (4.1) are the CHPT results given in Eqs. (|3.2| ) and (|3.4j) , 
whereas A p ^ ys denotes a phenomenological expression for the resonant amplitude: 

A% ys = A^pe tao . (4.2) 

The factor pe ia ° in the above equation takes into account the corrections to the tree- 
level amplitude of rj — > 37r°, which are known to be large [II]. Assuming a flat Dalitz 



Plot for this decay -not withstanding experimental constraints- and using the relation 
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B (m d - m u ) = m 2 K o - m 2 K+ - m 2 + m 2 + , (4.3) 

from the experimental data on T(i] — > 3tt°) [| we find p = 2.0 ±0.1. 

Contrary to p, the phase a® cannot be extracted from the rj — > 3tt° data. Similarly 
to the K — > 37r analysis of [Q, in order to evaluate ocq, we expand the one-loop 
rj — > 3tt° amplitude Elj around the center of the Dalitz Plot. Hence we obtain 



1/2 

«o = ^^(l-^| (2s + m 2 )^ 0.18, (4.4) 
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Figure 3: Partial decay rate of r] — > 7r°7r°77 as a function of the energy cut \slj^ — m^a \ < 8m. 
Full, dashed and dotted curves as in fig. 2. The two dash-dotted lines, denoting the absolute 
value of the interference between A p ^ ys and A]f^, have been obtained for a = 0.16 (upper 
line) and «o = 0.20 (lower line). 



where So = {jn^ + 3m^)/3. 

As anticipated, the two figures show clearly that in the non-resonant amplitude the 
one-loop result dominates over the tree-level one in the whole phase space. Moreover, 
for s 77 > 0.15m^ the non-resonant amplitude becomes non-negligible with respect to 
the resonant one. For s 77 > 0.20m^ the dominant contribution to the decay is provided 
by A NR . 

All the distributions have been obtained using p = 2 in A p ^ ys and considering the 
dominant tt-tt loops only. We have explicitly checked by means of Eq. ( |3.4| ) that the 
kaon loops give a very small contribution, confirming the statement of the previous 
Section that the 7r-7r loops dominate the non-resonant amplitude in the whole phase 
space. The dependence of our result on «o is quite small, as shown in fig. 3, while a 
precise determination of p is important, in order to determine with a good accuracy 
the region where the non-resonant amplitude dominates over the resonant one. Our 
analysis could be improved if high precision data on rj — > 3ir° were available. They 
would allow, in principle, to include the (small) D-wave contribution we neglected. 
However, even with the present uncertainty on the normalization factor p, we can 
state that the loop effect is in principle observable. Indeed, for s 77 > 0.23m^, we find 
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that the non-resonant contribution is larger than the resonant one by more than a 
factor of 2. 

Notice that our distributions differ from those in |l[ -apart from the one-loop 
corrections to A^r- by an overall normalization factor. Since we agree with the 
analytic expressions reported in [|TJ, this discrepancy can be traced to a problem in 
the program used to produce their plots. 



5 Discussion and outlook 

In this paper we have explicitly calculated the dominant one-loop corrections in 
CHPT to the decay r\ — > 7r°7r°77, going beyond the lowest order treatment of (ij 
- the latter essentially corresponds to a simple current algebra calculation. The 
phenomenological interest of this process is due to the experimental facilities which 
can effectively act, in the next few years, as 77-factories. A similar physical motivation 
led recently many authors to calculate leading corrections to the lowest-order CHPT 
prediction in both decays, such as rj —>■ 7r°77 [113, and scattering processes, such as 



77 — > 7r + 7T JT^| and 77 — > 7r 7r [15]. Earlier work on the CHPT predictions of 



pion polarizabilities can be found in [IB], and the comparison with forward-angle 



dispersion sum-rules is discussed in [fT7| . Recently the charged-pion polarizabilities 
have been computed to two loops flcfl . 

Recent results on 77 — > 7r°7r 7r° have inspired us, since there the lowest-order 
amplitude is suppressed and the corrections due to chiral loops dominate the cross- 
section. We found a similar result to hold for the non-resonant contribution to the 
decay r] — > 7r°7r°77. 

Despite the enhancement due to the one-loop corrections, the non-resonant ampli- 
tude is shadowed from the resonant (i.e. the 7r°-exchange) contribution, over a large 
portion of the diphoton spectrum. We have shown, however, that for large s 77 the 
one-loop corrections to the non-resonant amplitude dominate also over the resonant 
contribution (see figs. 2 and 3). A measurement of the partial width of 77 — > 7r 7r°77 
in this kinematical region -within the reach of the future facilities- would represent 
a new interesting test of CHPT at 0(p 6 ). 

We wish to conclude this Letter with a few comments concerning the possibility 
of future developments. In particular, one might ask why we concentrated on the 
neutral pion channel, rather than calculating the amplitude of i] — > 7r + 7r _ 77, which 
is statistically favored. The reason is that the decay rj — > 7r + 7r~77 is dominated by 
the bremsstrahlung of r] — > 7r + 7r~7 [|TJ. Since the latter is not suppressed already at 
the tree level, we expect that the one-loop corrections not related to the i] — > 7r + 7r~7 
amplitude will be hardly detectable. | From this point of view, when considering the 
outlook for a future calculation and a potentially related measurement, we expect 
that the scatterings 77 — > 7T + 7T~i] and 77 — > 7r 7r°?7 will provide more interesting tools 
for the study of chiral-loop effects. 



2 From the theoretical point of view, in order to isolate these effects it is necessary to implement an ap- 



propriate definition of 'generalized bremsstrahlung', similarly to what has been done in 1 10 for the radiative- 
four-meson amplitudes. 
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Added note 



After submitting this Letter, we learned about a similar calculation which included 
also the rj-rf mixing effect and an estimate of 0(p 8 ) contributions in the resonant 
amplitude ||19| . This result confirms that the effect we calculated is the dominant one 
in the non-resonant amplitude. 
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Appendix 

The general definition of C2o{x, y) in terms of the three-denominator one-loop scalar 
functions can be found in [|10j. In the n-n case and for x, x — 2y > Am\ the explicit 
expression is given by: 

(47r) 2 KeCVo(x,y) 



(16ir)$mC 20 (x,y) 



where 0o 
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